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Abstract 

The asymptotic scheme of post-Newtonian approximation defined for general relativity (GR) 
in the harmonic gauge by Futamase & Schutz (1983) is based on a family of initial data for the 
matter fields of a perfect fluid and for the initial metric, defining a family of weakly self-gravitating 
systems. We show that Weinberg's (1972) expansion of the metric and his general expansion 
of the energy- momentum tensor T, as well as his expanded equations for the gravitational field 
and his general form of the expanded dynamical equations, apply naturally to this family. Then, 
following the asymptotic scheme, we derive the explicit form of the expansion of T for a perfect 
fluid, and the expanded fluid-dynamical equations. (These differ from those written by Weinberg.) 
By integrating these equations in the domain occupied by a body, we obtain a general form of 
the translational equations of motion for a 1PN perfect-fluid system in GR. To put them into a 
tractable form, we use an asymptotic framework for the separation parameter rj, by defining a 
family of well-separated 1PN systems. We calculate all terms in the equations of motion up to the 
order r/ 3 included. To calculate the 1PN correction part, we assume that the Newtonian motion of 
each body is a rigid one, and that the family is quasi-spherical, in the sense that in all bodies the 
inertia tensor comes close to being spherical as r\ — > 0. Apart from corrections that cancel for exact 
spherical symmetry, there is in the final equations of motion one additional term, as compared 
with the Lorentz-Droste (Einstein-Infeld-Hoffmann) acceleration. This term depends on the spin 
of the body and on its internal structure. 
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I. INTRODUCTION 



Explicit, tractable equations of motion of celestial bodies can be used to compute 
ephemerides which have to represent the prediction of a given theory of gravitation for these 
motions. Such tractable, actually-used equations of motion are among the most important 
ones in a theory, of course. Naturally also, such equations are necessarily of an approximate 
nature, and this is all the more so in relativistic theories of gravitation, all of which are 
much more complex than is Newton's theory. Most of the parameters that enter these 
equations of motion are adjusted for the very construction of the ephemerides, or at least 
are adjusted under the assumption that the theory considered is correct . Therefore, 
the ephemerides represent in fact a fitting of astronomical observations by the given theory 
of gravitation or, more precisely, by the actually-used, approximate equations of motion 
derived from that theory. Einstein's general relativity (GR) is by very far the mostly 
investigated theory, and the literature about the post-Newtonian (PN) equations of motion 
for the mass centers (EMMC's) of an isolated system of celestial bodies in GR is quite vast. 
However, the EMMC's that are actually used in relativistic celestial mechanics 
were obtained two years only after the proposal of GR: these are the Lorentz-Droste (LD) 
equations j^J. They are more widely known as the Einstein- Infeld-Hoffmann equations, 



because that work of Lorentz and Droste 



4] 



seems to have been forgotten until rather 



recently [5|, and because these same equations were derived (though under different 
assumptions and for two bodies only) by Einstein and his coworkers 

The most satisfactory derivation of the LD equations is that provided by Damour, Soffel 
and Xu (hereafter "DSX" for short) [5|. As shown by DSX, the LD equations are the 
"monopole-truncated" EMMC's, corresponding to the ideal case where the gravitational 
field of each body is characterized by just one parameter (its constant mass), all higher-order 
multipole moments and all spin moments being zero. For each body, a set of such moments is 
defined by these authors generalizing results obtained by Blanchet and Damour (which 
were based on an exact multipole expansion of the retarded potential due to a compact 
source) Q], these moments are defined from multipole integrals, limited to this body, of 
the relevant current and energy densities; these moments then determine the coefficients in 
multipole series expansions of the (external) PN gravitational potentials produced by this 
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body. The results of DSX j^l Isj], in particular the exact recovering of the LD equations as 
the monopole-truncated EMMC's, are further confirmed by Racine and Flanagan 9] who 
also consider multipole expansions, but use a somewhat different approximation method 
(based on a separation between Newtonian and (first) post-Newtonian equations and, 
therefore, leading to Poisson equations instead of d'Alembert equations for the gravitational 
potentials; moreover, in order to extend the validity of the results to strong fields, Racine 
and Flanagan use surface-integrals definitions of the multipole moments and the mass 
centers). A question arises as to the relevance of the monopole-truncated model (or, shortly, 
the monopole model) leading to the LD equations: does this model approach the relevant 
exact general-relativistic motion accurately enough? (The fact that the LD equations do 
allow to compute accurate ephemerides does not give a strong argument for answering 
"yes:" from the purely logical point of view, one cannot preclude that the observational 
agreement could be less good if one would use a more accurate approximation.) Within 
the DSX framework the answer to this question should be obtained by incorporating 
higher-order multipole moments and spin moments, and by checking the magnitude of 
the corresponding new tams in th e e qUationS of m0ti on. DSX Q We indeed StU d,ed 
a "monopole-dipole" model for which, in addition to the monopole mass moment, also 
the spin dipole moments are non-zero, and in fact may be considered constant. Although 
DSX [8] did not discuss the magnitude of the new (spin-orbit and spin-spin) terms in the 
equations of motion, one may guess that they should be very small in the solar system, 
because they turn out to be of a high order in the separation parameter. However, in view 
of the infinite number of multipole and spin moments, there is no a priori guarantee that 
the set of the other moments gives a negligible contribution. Indeed DSX "state clearly that 
these models [such as the monopole model and the monopole-dipole model] do not need to 



represent the first step in some asymptotic approximation to reality, but only to be a_b 
the phenomena" with an acceptable accuracy, and in a logically consistent manner." 



e to "save 



By definition, neither in the monopole model nor in the monopole-dipole model does 
the internal structure of the gravitating bodies influence the equations of motion. However, 
there are indications that in fact the internal structure might play a non-negligible role. 
There is first a qualitative argument: in any relativistic theory, the mass-energy equivalence 
implies that any kind of energy should both contribute to the gravitational field and be 
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subjected to its influence. Thus, for instance, the rest-mass energy, but also the energy due 
to the interaction between matter and gravitational field, of which the Newtonian potential 
energy gives a first approximation, should play a role — as does the kinetic energy due to the 
internal motion: the influence of the latter is attested by the presence of the spins in the 
equations of motion derived from the monopole-dipole model. And since the distribution of 
the energy among these different forms can vary from one celestial body to the other, one a 
priori expects that the different distributions might affect the motion. A more quantitative 
argument follows from a derivation of the EMMC's in an alternative scalar theory, which 



has been done recently 
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12 1 . Indeed, it has been found there that several structure 



parameters enter the explicit EMMC's tailored to eliminate numerically-negligible terms, 
and it has been argued that the same should occur "in nearly any other theory" if a similar 
method was used 12] . 



Therefore, the goal of the present investigation was to check whether this conjecture 
about the influence of structure parameters does apply to GR. The method used in Refs. 

nnn 

[10, HH 1131 - which will also be used in the present work, is very different from that followed 
by DSXp, III and by Racine and Flanagan 9], who use multipole expansions from the 
beginning and who, in a first step, analyse the equations for the moments in reference 
frames attached to the different bodies. 1 Instead, our first step consists simply of an 
integration of the PN field equations for a perfect fluid in the global reference frame, say F, 
which provides a general (but not tractable) form of the PN EMMC's — the mass centers 
themselves being defined as local barycenters of the PN rest-mass density in the frame 
F Q. Such genera. PN EMMC's were not given in previons works, except in the form 
involving infinite series of multipole and spin moments 9], which is quite complex. Of 
course, our restriction to perfect fluids means some loss of generality as compared with the 
approaeh of Refs. ABB, whieh does not need to consider a partien.ar material model; 
but this restriction is not a serious one in the solar system 18], essentially because, inside 



1 Our method is closer to that used by Fock [l^, but it differs from it in that we use the "asymptotic" 
scheme of PN approximation \lA llfij , instead of the "standard" scheme developed by him [lj] and by 
Chandrasekhar 16]. The difference between these two schemes will be demonstrated in Sect. [H] especially 
Subsect. Ill CI In the particular case of a test particle in a Schwarzschild field, these two schemes are yet 
equivalent |l7| . Another important difference is that we use also a definite asymptotic framework for the 
separation parameter, which is exposed in Sect. IIVI 
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massive bodies, the deviatoric stresses are small as compared with the isotropic pressure. 
(Moreover, the present method might be used with a more sophisticated material model, 
at the price of redoing the calculations.) In a second step, to get tractable equations, we 
use three simplifications which do occur for the solar system: i) the fact that its main 
bodies are well-separated. This is an essential assumption, which is set by assuming that 



a precise separation parameter [111 ], which we denote by 77, is a small number 770 for the 
system of interest, S, and by introducing a family (S v ) of well-separated systems (each 
system is weakly gravitating, with the field strength being nearly the same for all systems, 
and with S 770 = S). It enables one to define a hierarchy of well-defined approximations for 
the "tidal" effects, by calculating asymptotic expansions, with respect to rj, of the integrals 
that enter the general form of the PN EMMC's jljj]. ii) The second simplifying feature 
of the solar system is that its main bodies have a nearly rigid rotation about themselves, 
iii) Lastly, the main bodies of the solar system are nearly spherical. We use the sphericity 
assumption only for the zero-order rest-mass density, po (but this implies that the zero-order 
pressure and the self Newtonian potential are also spherically symmetric in each body). 
But, as well as for the rigidity assumption, we use the sphericity assumption merely at 
the stage of calculating the PN corrections, not for the zero-order calculations themselves. 
Moreover, we shall simultaneously derive an explicit (though less condensed) form of the 
equations of motion without using this sphericity assumption, and using instead a much 
weaker "quasi-sphericity" assumption, expressed by Eq. (j4.11|) . 2 We argued previously 
[ill ] that, in the solar system, the sphericity assumption is likely to lead to relative errors 
smaller than 10~ 3 when calculating the PN corrections, which are already very small. In 
summary, we believe that the assumptions introduced in the present paragraph are enough 
justified in the solar system. The "asymptotic" PN scheme, which is used here, leads to a 
separation between zero-order (Newtonian) and first-order (1PN) equations. Therefore, the 
small effects neglected, such as the departure from rigidity or the shear stresses [3], could 



2 Thus, the incompatibility between exact spherical symmetry and non-zero self-rotation is not relevant here, 
firstly because we do give the formulas for the non-spherical case. Secondly, when we use the assumption of 
spherical zero-order densities to get a more tractable expression of the PN corrections, this assumption has 
the status of a relevant approximation. According to this view, the zero-order calculations might take into 
account the departure from sphericity. The spherical zero-order density po considered when calculating 
the PN corrections might then be defined as some relevant approximation (a least-squares approximation, 
say) to the exact, non-spherical zero-order density. The same can be said about the rigidity assumption. 
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be described accurately enough by taking them into account merely for the Newtonian 
calculations. (Such Newtonian calculations using a more general model of the bodies are, 
of course, a well-developed subject in geophysics.) Hence, due to this separation in the 
scheme used here, it would not make much practical sense to include such small effects into 
the calculations of the PN corrections, unless at the same time one would compute second 
post-Newtonian corrections to the motion. 

This paper is organized thus: the next Section summarizes the asymptotic-expansion 
method of the local fields and equations for a perfect fluid in GR under the harmonic 
gauge. The general form of the equations of motion for a perfect-fluid system in GR (in the 
harmonic gauge) is derived in Sect. 11111 The framework which is used to rigorously account 
for the good separation between bodies is summarized in Sect. IIVI Based on Appendix 1X1 
in which we compute the integrals entering the general form of the EMMC's, for the case 
of well-separated, rigidly-rotating and quasi-spherical bodies, Sect. fVl presents the explicit 
EMMC's obtained for that case, together with their specialization to the case with spherical 
mass densities. Our conclusion is presented in Sect. IVII 

II. ASYMPTOTIC P.N. APPROXIMATION FOR A PERFECT FLUID IN G.R. 
UNDER THE HARMONIC GAUGE 



The standard PN approximation is based on the classical works of Fock Il3| and Chan- 
. At those times, the numerical analysis of partial differential equations 
(PDE's) was far less developed than it is now. It would now seem relevant to build an 
approximation scheme in agreement with the asymptotic schemes currently used in the nu- 
merical analysis of PDE's. In a such scheme, one has a regular family of systems, (S e ) (each 
system being defined by a boundary-value problem — an initial-value problem for that mat- 
ter, because relativistic gravitational equations are hyperbolic), so that one may make the 
"small parameter" e as small as desired, indeed. Then, for the corresponding family of fields, 
one should be able to state some asymptotic expansions with respect to e; in particular, all 
of the unknown fields should be expanded, because in general all depend on e. Moreover, 
the family should be associated in a physically natural way with the system of interest S, 
which should correspond to a small value eo of the parameter, S = S eo , so that it makes sense 
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to use the expansions for the given system S. Finally, e should be the natural field-strength 



parameter 



13j . so that a small value of e means indeed a weak gravitational field. It has 



been shown by the work of Futamase and Schutz Q| how to develop a post-Newtonian 
approximation along this line for GR (in the harmonic gauge). Below, we indicate first 
(Subsects. Ill Al and 111 Bj) how, starting from their initial condition, one indeed obtains the 
general expansions and expanded equations derived by Weinberg in Sections 9.1 and 9.3 of 
his classical book j^J. Our aim is not to present mathematically rigorous proofs but to show 



that the main features of Weinberg's approach occur rather natural 



y and convincingly from 



the point of view of asymptotic analysis. Futamase and Schutz [14] did not investigate the 
relation between their "asymptotic" PN approximation and Weinberg's expansions, which, 
in our opinion, remained difficult to understand. Then we show (Subsect. Ill Cj) that, how- 
ever, the explicit expansion of the energy-momentum tensor T and the explicit expanded 
equations obtained according to this "asymptotic scheme" for a perfect fluid are different 
from the equations stated by Weinberg in his Section 9.8. 

A. Expansion of the matter fields and the metric 

Futamase and Schutz assume the following initial data for the fields p (pressure), /i* 
(proper energy density in mass density units), 3 u (coordinate velocity), and g [space-time 
metric; in this paper, we follow the conventions of Weinberg [2Q], in particular the signature 
is (3,1)]: at time t = 0, 

p (e) (x) = eV 1} (x), (2.1) 
/i^(x) = e y«(x), (2.2) 
u (£) (x) = £u (1) (x), (2.3) 



g® ij (x) = $iv (v /Z ^^ )iJ ) Q (x)=0 (l<i<3, l<j<3) (2.4) 



3 Futamase and Schutz |l4|. as well as Weinberg [20|, use the notation p for /z*. Moreover, they set c = 1, 
so that the physical dimensions are not apparent in their equations. 
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[here g = det(g^), (g^ u ) is the inverse matrix of (g^u), and x = (x l ) is the spatial po- 
sition]. Conditions (j2.1|) - ()2.3j) for the matter fields are suggested by an exact similarity 
transformation valid for the Euler-Newton equations, which defines the weak-field limit for 



a perfect-fluid system in Newton's theory jl4L Il5l|. This indeed suggests to define the initial 
conditions for the Newtonian limit of any "relativistic" theory by applying the similarity 
transformation to the initial data, at least for the matter fields — and this leads exactly to 
Eqs. (j2.1|) - ()2.3j) . As to condition ()2.4|) . it is just the condition that one would impose if one 
would wish to have the following "conformally-Euclidean" (or isotropic) form for the space 
metric 7 (the spatial part of g in the frame defined by the coordinate system): 

Hi = (2-5) 

without imposing any a priori restriction on the factor yj—g. (Hence, it is legitimate to 
postulate this condition (|2.4jl . but it is not clear that one has to postulate it, unless one 
does so precisely to enforce spatial isotropy.) Note that condition ()2.4|) does not depend 
on e. Thus, using first Eqs. (|2.1|) - (j2.3j) to go from the finite small value e , valid for the 
physically given system S, to e = 1, and then once more to go from e = 1 to the arbitrary 
value e, one indeed naturally associates with S a family (S e ) of systems. 

Now, let us change the mass and time units for system S e in this way jl^: [M] e = e 2 [M] 
and [T] e = [T]/e, where [M] and [T] are the units for system Si. In these units, the initial 
data (j2.4j) for the metric is unchanged (because the metric is adimensional) , but Eqs. (j2.1j) - 
(12. 3|) become simply 

p {e) (x) = (x) , ^ (x) = (x) , u (e) (x) = u (1) (x) . (2.6) 

Thus, the initial condition is independent of e, moreover we have just e = c _1 (we take the 
velocity of light to be c = 1 in the units for system Si). The initial- value problem still 
depends on e or more exactly on e 2 , since in general units the Einstein equations involve the 



square c 2 = e 2 : 



Rtw _ \ g ^ R \ = -^—T" V (2.7) 
2 c 2 



(this writing is in accordance with Weinberg's, which itself follows his conventions about the 
Riemann and Ricci tensors; we take T in units of mass density, thus ML -3 ; note that G with 
dimension L 3 M _1 T -2 is invariant in the change of units). The effective small parameter is 



thus e 2 , rather than e. Hence, it suggests itself to state Taylor expansions with respect to 
e 2 = c -2 , starting from the zero-order term (denoting with a prime the fields as expressed 
in these varying units): 



P' = Po+Pi c +0{c 



fjf* = ^ * +yU ' 1 * c -2 + o( c - 4 ), u ' = u' + u' lC - 2 +O(c- 4 ), (2.8) 



etc. When coming back to the fixed units (of system Si), these expansions are modified in 
a straightforward way: 

p=[p +Pie 2 + O(e 4 )]e\ /i* = [^ + ^ e 2 + 0(e 4 )]e 2 , u = [u + u l£ 2 + 0(e 4 )]e. (2.9) 
On the other hand, the definition of tensor T from the matter fields is for a perfect fluid 



= [ji* +pc~ 2 )UW +pc- 2 g^ 



(2.10) 



13j,|20j], where is the four- velocity of the fluid: 



U" = dx^/ds, ds = (-g^dafdx") 



(2.11) 



hence 



IP = U°u l /c 



(2.12) 



since u % = dx l /dt and x° = ct. Therefore, the definition of T from the matter fields does 
involve directly c _1 = e. However, the equation Tf^ = 0, once written in terms of the 
matter fields using this definition, will contain only c™ 2 . (This is not immediate to check on 



the usual forms of the exact perfect-fluid equations in GR, e.g. Font 



it is set c = 1; see Ref. 



2l| . because usually 



lq . Eqs. (3.8) and (3.10), for the equations in an alternative 



theory, based on a slightly different conservation equation.) Thus, this passage by the 
HHHdnUnB. ^ least suggests that the expats in the usual fie.d-streugth parameter 
231 ] . or the formal expansions in 1/c |16(, could be in fact expansions in 



13, 



20, 



22 



the square A = e 2 , or 1/c 2 . (One sometimes invokes the invariance by time-reversal, but 
this applies to the equations, not to their solutions for a generic situation.) As stated by 
Weinberg j^j, "the real justification for these expansions will come below when we show 
that they lead to a consistent solution of the Einstein equations." We will use below the 
varying units for the metric, whose expansion is less straightforward. 
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In order that there are well-defined zero-order fields and equations, it is necessary that 
the metric has a finite limit as e — > 0. This should occur in the fixed units, for we know that 
in a real weak field, thus for a finite small value €q of the field strength, we can approximate 
the metric by a finite metric which is independent of eo, namely a flat one. Thus, without 
anticipating the flatness, the metric should have an expansion with a first term of order zero: 

g = g + O(e 2 ). (2.13) 

Setting 77 = diag(— 1, 1, 1, 1) and h = rj — y/—gg, the harmonic gauge condition writes 
h^f = and in this gauge the Einstein equations take the ("relaxed") form Q| 

Uh aP = -ff v h^ v = 16vrA Q/3 , (2.14) 

where A a/3 is the sum of a term linear in tensor T and a polynomial P af3 (h fiu , h^, h^), 
having at least quadratic terms. The zero-order equations are hence simply 

Dhf = lQnAf (2.15) 

(with and Aq^ the zero-order coefficients in the expansions of h alS and A af3 ). From 
the expansions ([2.9)1 . it results that tensor T is 0(e 2 ) (in the fixed units) and hence makes 
no contribution to A^. Hence, we have A^ = P a/3 (ho U , h^ p , ho"). Thus, the zero-order 
equations ([2.15)1 are still nonlinear, just as the equations for the exact field h. However, 
ho = (hence go = rf) is a solution of these equations which is compatible with the initial 
data deduced from (|2.4|) . h(t = 0, x) = and doh(t = 0, x) = 0. Moreover the solution 
of the initial-value problem for these nonlinear wave equations should be unique. Thus the 
limiting metric is flat indeed: 

go = T). (2.16) 

However, if one reexpresses this limit in our varying units (or, for that matter, if one just 
uses the varying unit of time, thus replacing the time t by the "dynamical time" t' = et 
ji3 J3) he finds that the limit go has the following components, one of which is singular at 
the limit e — > 0: 

(o^J = diag(-e- 2 , 1, 1, 1) = diag(-c 2 , 1, 1, 1). (2.17) 

From 1)2.17)1 . it follows that, to 1PN order, i.e., including terms up to the order c -2 in the 
varying units, it is natural to postulate the following expansions: 

9oo = -c 2 + iSoo + 2^00 c- 2 + 0(c" 4 ), (2.18) 
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9' ij ^S ij+1 g' ij c- 2 + 0(c- 4 ), (2.19) 

SL= i^c- 2 + 0(c- 4 ). (2.20) 
This leads, in the fixed starting units, to 

3oo = -1+ i£oo e 2 + 2^00 e 4 +0(e 6 ), gij = <%+ igij e 2 +0(e 4 ), g 0i = ig 0i e 3 +0(e 5 ). (2.21) 



This is the expansion written by Weinberg |2£j] (accounting for the fact that this author 
incorporates the small parameter in the expansion coefficients) and used by Racine and 
Flanagan p. 

The asymptotic expansions of the fields have first of all to be valid at fixed values of 
the time and space variables, of course. Since the velocities in the system S e vary like e, it 
follows that the characteristic times (e.g. the orbital periods) vary like e~ l . Therefore, the 
relevant time variable in the expansions is not the time t in fixed units, but the "dynamical 



time" t' = et 



ay 



151 ]. This is in fact obvious if one remembers that the initial data in the 
varying units (and thus with the time t') is independent of e. Thus, the coefficients of the 
expansions: po, p±, etc., are functions of x and t' but they are by definition independent of 
e. With this in mind, the insertion of ()2.9j) into the initial data (|2.1jl - (|2.Hjl gives the initial 
data for the coefficients: at the initial time, 



p (x) = p 



«(x), pi(x)=0, (2.22) 



^(x)=// (1) (x), /^(x)=0, (2.23) 



u (x) = u«(x), ui(x) = 0. (2.24) 
B. Expanded field equations 

It is easy to check that the calculations done by Weinberg in his Sections 9.1 and 9.3 
follow exactly from the expansions ()2.9j) and (|2.21j) of the matter fields and the metric. (We 
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stay in fixed units with c = 1 until the end of Section |TTJ) One difference is that, because 
the relevant time variable in the asymptotic expansions is the dynamical time if = et, the 
coefficients in the expansions, e.g. igfy in Eq. f|2 . 2 lj) o . are functions of t' (and of the space 
coordinates x k ). Hence, when differentiating the e-dependent fields with respect to x° = t, 
as this occurs e.g. in the definition of the connection and the Ricci tensor, one increases 
of one the order in e (which is what Weinberg essentially assumes), but also one ends up 
with a derivative of the expansion coefficient with respect to t' . For instance, in uniform 
conditions, differentiating Eq. (|2.21Jl o gives 

% = %^ 3 +°( £5 )- (2 . 25) 

Another difference is that the unique small parameter e now appears explicitly in the ex- 
pansions (but of course it does not in the expanded equations, which obtain by coefficient 
identification in the expansions). Thus, for instance, the components Rqo, Rio and o f 
the Ricci tensor are indeed of order 2, 3 and 2 in e respectively, Eqs. (9.1.23-25) in Ref. [20], 
but (at least for a fluid), the corresponding components T 00 , T*° and T u ' of tensor T are 
now respectively of order 2, 3 and 4, in the same e, hence the expansions 

= e n »» ( T^ + X T^ e 2 + 0(e 4 )) , n Q0 = 2, n i0 = 3, n l3 = 4. (2.26) 



{Compare Eqs. (9.1.42-44) of Weinberg [20j .} Defining tensor = — ^g^T^, one then 
finds that 5*00, S i0 and S^- do have expansions with the same orders as R 00 , R i0 and ir- 
respectively, the coefficients of these expansions being given by his Eqs. (9.1.49-52). The 
Einstein equations (J22J), equivalent to R^ u = SttGS^, thus split to Weinberg's equations 
(9.1.53-56). It follows that the expansion (|2.21j) of the metric is explicitly [his Eqs. (9.1.57- 
58,60-61,63-64)]: 

g 00 = -l- 2$e 2 - 2($ 2 + ^ + 0(e 6 ), 9ij = %(1 - 2$e 2 ) + 0(e 4 ), 9oi = 2Q e 3 + 0(e 5 ), 

(2.27) 

the potentials $, ip and Q being solution of the Poisson equations 

A$ = 4ttG T 00 , (2.28) 



02 ,T, 

A ^ = df2 + 4 ^(iT 00 + T il ), (2.29) 
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AQ = lQnG T i0 , (2.30) 

where A is the Laplacian corresponding to the Euclidean metric which has components 5ij 
in the given harmonic coordinate system utilized, thus Aip = ip. 



M ■ 



In the same way, starting from the general expansion ()2.26j) of the energy-momentum 
tensor T and the explicit expansion ()2.27j) of the metric g, one checks straightforwardly 
that the same expanded equations are derived from the local dynamical equations TJ£ V = 
as Eqs. (9.3.2-5) in Ref. j20|, namely 

d( T^) a( T^) 0* 

~m^ + ^x^-' oT m' (2 - 32) 

d{iT o 0) d{iT o 3) * $ 

+ ~ oT w> (2 - 33) 



^(i^ 0t ) , d(iT^) _ 00 9$ _ 00 
dt' dxi 1 dx i 



g(2$ 2 + ^) | dd 
dx i dt' 



+ oT o j{ dCi_dC i 



dx l dxi 



n- 9$ 9$ 9$ 

Note that, so far, the assumption of a perfect fluid was necessary only to ensure that the 
components of tensor T have orders in e given by Eq. ([2.26)1 . which can be easily checked to 
be consistent with the expansions ()2.9j) of the matter fields and with the definition of tensor 
T for a perfect fluid, Eq. (|2.1(J|) above. (The next Subsection shows the explicit expansions 
of tensor T that one gets thus.) Of course, much more general constitutive laws would 
still lead to the same orders for tensor T. One would then have to adapt the initial data 
(|2.1|) - ([2.3p and the expansions ()2.9|) correspondingly, however. 



C. Explicit expansion of the energy-momentum tensor and expanded equations 
for a perfect fluid 

To obtain the PN field equations for a perfect fluid, we must first write the explicit 
expansion of the energy-momentum tensor of a perfect fluid as function of the expansions 
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of the matter fields and the metric field. We insert the matter fields expansions ()2.9|) and 
the expansions of g^ u and U° deduced from (j2.21|) by Weinberg, into the expression of the 
energy-momentum tensor for a perfect fluid (|2.1U|) . in which c = 1. This leads easily to the 
explicit expansion of the energy-momentum tensor: 

T 00 = ^, 1 T 00 = ^+/iSK-2$), (2.35) 



T Oi = ^u l , xT ' = K + M - 2$ ) + PoH + (2-36) 



(2.37) 

This expansion is different from that obtained by Weinberg {Eqs. (9.8.4-6), (9.8.10-13) 
in Ref. |20j}. This is because Weinberg does not expand the matter fields of a perfect 
fluid: p, jj* and u. This same difference between the standard PN scheme of Fock [.13J and 
Chandrasekhar 



lq and the asymptotic expansion method has previously been noted by 
Futamase and Schutz ^| and by Rendall [24J]. Weinberg's approach is distinct from those 
of Fock and Chandrasekhar in that, in a first step, Weinberg derives the PN equations for a 
general form of the energy-momentum tensor, which tensor he does expand. For this reason, 
Weinberg's equations essentially coincide with those got with the asymptotic scheme (the 
differences being explained at the beginning of Subsect. Ill BJ) . until the explicit expansion 
of tensor T for a fluid is written. From this stage, Weinberg's equations [his Sect. (9.8)] 
coincide with those of Chandrasekhar 16], which are essentially equivalent to those of Fock 
j"!^ . In Weinberg's work, the departure from the asymptotic scheme can be seen e.g. in the 
fact that the same field denoted p by Weinberg cannot be interpreted at the same time as the 
exact proper energy density, which we note /x*, as it is in his Eq. (9.8.1), and as the lowest- 
order coefficient in the expansion of T 00 , thus our /Xq, as it is in his Eqs. (9.8.4) and (9.8.9). 
Moreover, Eq. (|2.35j) 9 here, that gives the following coefficient in the expansion of T 00 in the 
asymptotic scheme, has the correction term * as compared with the c—ding equation 
of Weinberg, his Eq. (9.8.10). As suggested by Futamase and Schutz (14J, the interpretation 
of the unexpanded matter fields of the standard scheme that is closest to bridging the gap 
with the asymptotic scheme, is to consider them as the second approximations of the exact 
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fields. This interpretation means assuming the following correspondence between Weinberg's 
notation in his Eqs. (9.8.4-15), and the present notation: 

P A*(i) = e 2 (/^o + V*i e2 )i P ^P(i) = e4 (Po + Pi e2 )> v ^ u (1) = e(u + uie 2 ). (2.38) 

With this interpretation, his Newtonian equations [Eqs. (9.8.4-9)] are valid up to 0(e 2 ) terms 
not included, and his 1PN equations [Eqs. (9.8.10-15)] should be valid up to 0(e 4 ) terms 
not included — provided that, when needed, one inserts the relevant power of e, the small 
parameter being incorporated in the coefficients in Weinberg's notation. Thus, for instance, 
Eq. (9.8.4) of Weinberg becomes 

r 00 e 2 = A iJ l) (l + O(c a )) > (2.39) 

and his Poisson equation (9.8.9) for the Newtonian potential $ becomes 

e 2 A$ = 47rG/i^(l + 0(e 2 )), or A$ = 4nG(n* + /^e 2 ) + 0(e 2 ), (2.40) 

which is obviously compatible with the exact equation of the asymptotic scheme, Eq. ()2.28|) 
with (f2~3H|b: 

A$ = 47rG/i*. (2.41) 

However, if one uses Eq. (|2.40|) instead of the exact equation (|2.41|) to compute then 
$ is determined only up to unknown 0(e 2 ) terms, really. (Indeed, except at the initial 
time, Eq. (|2.23jh the coefficient /i* will of course not be zero.) But the Newtonian potential 
<3> intervenes already at the lowest order in the 1PN equation of motion [Eq. (9.8.15)], 
namely by the term — pV$ in this equation. Hence, if one does interpret Weinberg's 
unexpanded matter fields as the second approximations of the exact fields, thus writing 
Weinberg's equation (9.8.9) as Eq. (J2.40)) above, then his 1PN equation of motion (9.8.15) 
is accurate only up to unknown 0(e 2 ) terms, in fact — i.e., it is not more accurate than the 
Newtonian equation of motion. This is the reason why, in our opinion, the standard PN 
scheme is not compatible with the asymptotic scheme. We note that the scheme used by 
Damour et al. for 1PN approximation p, ls| does not pertain to the standard scheme of 
Fock and Chandrasekhar, because, in the DSX scheme, neither the gravitational field nor 
the matter fields are expanded: DSX consider the second-approximation fields, which we 
note U(i),p(i), (i)^, etc., and they do not split them into zero-order and first-order parts. 
In the standard scheme, the gravitational field is expanded (split), but the matter fields 
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are not. In the asymptotic scheme, all fields are expanded. Hence, we feel that the DSX 
scheme is compatible with the asymptotic scheme, though it differs from the latter. As to 
Racine and Flanagan they use that part of Weinberg's approach {§§ (9.1) and (9.3) in 



Ref. 



20j]} which is compatible with the asymptotic scheme. 



Our aim in the present paper is to follow the asymptotic scheme until tractable equations 
of motion are got. To do that, we will moreover assume barotropic perfect fluids (one fluid per 
astronomical body). Recall first that the proper energy density p* that enters the expression 
(I2.1()j) of tensor T for a perfect fluid is, precisely, the sum of the proper volume densities of 
rest-mass and of (elastic) internal energy, p* = p*(l + II). The assumption of a barotropic 
fluid means that p* = F(p) depend only on the pressure p, as well as does II, the latter 
being given by |l3j] 

For a barotropic fluid, it is hence convenient to replace the initial condition ()2.2|) for the 
proper energy density p* by one for the proper rest-mass density p*\ 

at t = 0, p* (e) (x) = e 2 p* (1) (x) instead of /i* (e) (*) = e 2 /i* (1) (x), (2.43) 

which makes p* order 2 in e, like p*. Accordingly, the initial condition for the expansion 
coefficients f!2.23|) is replaced by 

att = 0, p*(x)=p*W(x), pj(x) = 0. (2.44) 

iously Q] 

HQ 



To ensure that the pressure and density fields obey Eqs. (|2.1|) and (|2.43|) simultaneously, one 
assumes that the function defining the barotropic state equation for system S e is 



F e (p) =e 2 F 1 (e- 4 p). (2.45) 

The expansions of the auxiliary matter fields p*, II, and p* follow then, of course, from the 
expansion (|2.9jb of the leading field p jisl ]: 

p* = e 2 ( P ; + p* 1 e 2 + 0(e A )), f^ = F 1 {p ), p\=p x F[{p Q ), (2.46) 



n = e 2 (n + O(e 2 )), U = G 1 (p ), (2.47) 
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p* = p* , n^pZUo + p*. (2.48) 
Since Eq. f!2.42|) is equivalent to the isentropy equation: dll = — pd(l/p*), a barotropic fluid 



is isentropic; hence, as shown by Chandrasekhar 



251 ]. its rest-mass is exactly conserved. This 



suggests introducing the density of rest-mass with respect to the Euclidean volume measure 



dV = dx dx dx associated with the considered harmonic coordinate system 



„* 



(2.49) 



which obeys thus the usual continuity equation with the velocity u. The expansion of p is 



p = e 2 (po + pie 2 + 0(e 4 )), p = p* = p*, Pl = p\ + p (-3$ + u 2 /2). 



(2.50) 



Using (ESDI) 2 in P^5]i. %2~3^ i and (|2~H7|li. and reporting in Eqs. (|2~3l>(|2~32]) . one checks 
that the latter ones reduce to the Newtonian equations in which p plays the role of the 
Newtonian density: 

dep Q + d j (j>o4)=0, (2.51) 



dfiPoul) + djipoUQU^) = -p $ i - Po,i- (2.52) 

Thus, Eq. f|2.31|) {Eq. (9.3.2) in Ref. |20(} expresses the zero-order conservation of mass, 
not the PN mass conservation. But, rewriting the {T^ 's [Eqs. f!2.35|) 0- 1)2.37^ ] in terms of 
Po and pi, one gets ()2.33|) as 



d t >(w + pi) + dj[(w + po + pi)u 3 Q + p u{] = po dt>$, w = p ( y + $ + n J , (2.53) 

which, combined with the Newtonian energy equation deduced in a standard way from ()2.51|) 
and ()2.52|) . gives 15|,|24j the order-one component of the continuity equation: 



d t 'Pi + djipiul + p Q u{) = 0. 
It remains to find the field equation for 1PN correction to the fluid motion. Setting 



#i = pi + po (^-3$ + n ^ +p , 



(2.54) 



(2.55) 



o x = xT 00 + T jj = pi + po ( + $ + n ) + 3p , 



(2.56) 



18 



we may first rewrite the expansion coefficients ()2.36|) 9 and (|2.37|) 9 as 

xT°* = Po u\ + 9x4 + 4( T 0i )$, (2.57) 

xT ij = p 4u{ + p u\ui + 9x44 + (p! - 2p <S> )6 ij + 4( T^>. (2.58) 

Inserting this into Eq. ()2.34|) . and simplifying terms with Eq. ()2.32|) . we obtain the sought-for 
equation: 

d t '{pau\ + 9x4) + dj(po4 u i + Po44 + °i44) + d i{Pi ~ 2 Po$) 

= -ax$,i - p [iJ,i + fyCi + (d,k ~ Ck,i)4}- (2-59) 

III. GENERAL FORM OF THE 1PN EQUATIONS OF MOTION FOR A 
PERFECT-FLUID SYSTEM IN G.R. (IN THE HARMONIC GAUGE) 

A. The definition of the mass centers and its motivation 

In GR, as in any relativistic theory of gravitation, any form of material energy must 
both contribute to the gravitational field and be subjected to its action. It is not obvious, 
therefore, to state which energy density may be used as a weight function so as to define 
relevant mass centers. Two arguments justify the choice made [l0[, of the rest-mass 
density p in the global reference frame, Eq. ()2.49|) : i) since this density obeys exactly 
the usual continuity equation, one may commute time differentiation and barycentration. 
(In other words, the velocity of the mass center of the body equals the average velocity in 
the body.) This was advocated by Will practical advantage, which it certainly 

is. In our opinion, this property is also an essential feature of the classical (Newtonian) 
definition of the mass center, without which the notion of a mass center becomes less useful 
physically; for the loss of this property means that, following the motion of all constituents 
of the body, one could not tell what is the motion of the mass center, ii) Rest-mass 
is well-correlated with astronomical observations, because it is indeed the presence of 
matter in the usual sense, thus characterized by its rest-mass density, that leads to the 
electromagnetic emission detected by the telescope. In practice, one determines an "optical 
center" which, once corrected from the "phase effects," is used to define the astronomical 

nn 

(observational) position of the mass center |26l . |27| . Hence, when giving the theoretical 
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definition of the latter, one has indeed to check that the density selected is well-correlated 
with the luminous density. Thus, in our opinion, the mass centers defined with the help 
of p are directly relevant to the astronomical observations. Several works have also used 
the choice of the rest-mass density (|2.49j) as the weight function to define the mass centers, 
e.g. Fock 0, Will Q, Brumberg |28[. (These works were based on the standard PNA 
scheme, however, which is distinct from the asymptotic scheme used in the present work: 
see Sect. In the literature on the PN equations of motion in GR, it has been more 

usual to define the mass centers as the local barycenters of a density p' which is obtained 
by adding three (small) densities to p: the density of internal energy, the energy density 
associated with the self Newtonian potential of the relevant body, and the density of the 
kinetic energy associated with its motion with respect to a local frame attached to the mass 
center of that body. See e.g. Misner et al. 22|, Spyrou H, Will H {cf. Eqs. (6.21) and 



(6.25) in the latter work}. Instead of the continuity equation, that density p' would obey 
a balance equation with a source term (this equation is not written in the quoted works 
[22! 0, 13 )■ m the DSX formalism, the definition of the mass centers is more involved 
{see around Eq. (5.10) in Ref. [3]}, but it also is not based on a density obeying the usual 
continuity equation. For sure, i) does not hold true with such choices, although the final 
equations of motion may be simpler. In addition, such densities as p' might be slightly 
less well correlated with the luminous density, because there is no reason that the latter 
increases with, for instance, the local density of the proper kinetic energy. However, the 
three additional densities are, of course, very small [0(c~ 2 )] as compared with p. More 
precisely, they are of the order p x GM' a j (c 2 r a ) where M' a is the total mass-energy in body 
(a) and r a is its size. Therefore, by itself, the difference in the densities p and p' which may 
be chosen as the weight function can imply only small and non-secular differences in the 
positions of the mass centers, at least in a weakly-gravitating system such as the solar system. 

According to DSX 5], one should consider carefully-chosen local reference frames to 
properly analyse the "internal problem" of celestial mechanics (i.e., that of determining the 
motion of each body around its own mass center), in order to avoid introducing metrical 
effects due to the global velocity and the external potential. (See also Kopeikin and Vlasov 
(3^].) In the present work, we limit the PN calculations to those necessary to get the 
translational equations, thus to solve the "external problem" in the terminology of DSX. 
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Using the expanded field equations of the asymptotic scheme, derived in the foregoing 
Section, we shall derive in this Section a general form of the 1PN translational equations 
for the "p-centers," Eq. ()3.17j) . The remaining work to solve the external problem at the 
1PN level is just to compute the integrals entering that equation. We shall compute these 
integrals up to and including the terms of the order 3 in the separation parameter (see 
Sect. II VI for a discussion of the good separation and the way we account for it). By using 
the simplifying assumption of a rigid Newtonian motion for each body [Eq. (|4.5[) ]. the 
only dynamical equations that we still have to use for that purpose are the Newtonian 
equations [in addition to the field equations for the gravitational potentials, Eqs. (|2.28j) . 
(12.29)) and (|2.3U|) ]. The Newtonian dynamical equations, when they are used, are written in 
the global reference frame: we do not have to introduce any "body-attached" reference frame. 



Thus, we define the exact masses and mass centers through the rest-mass density p 
[Eq. (jZI31)]: 



M a = / pdV, M a a EE / pxdV (3.1) 

where D a is the (time-dependent) domain made of the spatial positions x ee (x*) of the 
particles constituting body (a) (a = 1, ...,N) in the considered harmonic coordinate system 
(x M ). At the 1PN approximation, the mass and the mass center are approximated by 4 

M« =M° a +M l Jc\ M° a = f Po dV, Ml= f Pl dV, (3.2) 
M«a (1 ) ee I P(1)X dV = M°a + M^/c 2 , (3.3) 



with 



M a °a ee f Po xd^, Mlm ee f Pl xdV. 
JD a Jn a 



(3.4) 



Note that M° and an are the Newtonian mass and mass center. Using Eqs. ()2.51|) and (|2.54|) . 
one shows easily jlflj: i) that M a ° and Ml are constant in time (exactly so, insofar as the 



4 Henceforth, we shall use the expansions written with the "varying units" of mass and time, [M] e = e 2 [M] 
and [T] e = [T]/e, hence the expansion parameter is e 2 = c~ 2 . Note that the expansion coefficients and 
the expanded equations derived in Sect. [H]are valid in any units (although the expansions themselves 
are somewhat different in fixed units and in the varying units). What we are really interested in, is 
the physically given system S, corresponding to the finite small value eo of e. Therefore, the distinction 
between the times t and t' — et has no interest any more and we shall omit the prime. 
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matter fields cancel on the boundaries <9D a — in fact they are negligible there), and ii) that 



a d (l) 



/ p(i)U(i)dV + 0(c 4 ), u ( i) = u + uic 2 , (3.5) 

which is, at the 1PN level, the commuting property of time differentiation and barycentra- 
tion, referred to at the beginning of this Section; in direct connection with this, we have also 



/ ( Pl u + p ui)dy = M a 1 a 1 + 0(c- 2 ) (3.6) 

(this will be used later). Moreover, one finds from (|3.2)l - (|3.4j) that the PN correction to the 
position of the mass center is given by 



M 1 

a (D " a o = -YTTo ( a i " a o) + °(^ 4 )- ( 3 - 7 ) 



In the final equations of motion, we shall also use the notation 

x 0a = a 0; x la = C 2 ( a (l) — a o); x a = a (l) = x 0a + x laC 2 , (3-8) 

v a = a (1) = x a = x 0a + xi a c" 2 . (3.9) 
B. General form of the 1PN equations of motion 

As mentioned in the Introduction, the (1PN) equations of motion of the mass centers 
(EMMC's) are got by integrating the spatial components of the 1PN field equations of 
motion in the domains D a . Since the equations for the orders zero and one are separated 
[Eqs. (|2.52|) and (|2.59|) ]. the same occurs for the EMMC's. The integration of the zero-order 
field equation gives simply the Newtonian equation of motion 

K°a* =- / po^fdV, (3.10) 



where an upper dot indicates time derivative, and where we use Fock's 13] decomposition 

$ = $ (a) +0 a , (3.11) 

$<«)( X ) = _ G r f Po (y)dV(y)/|x-y|, (3.12) 

a (x) = —G f p (y)dny)/|x- y | (3.13) 
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(the corresponding usual expression of <3> being the solution of Eq. (J2.28JI under the usual 
boundary condition). When integrating the field equation ()2.59|) for the PN correction, we 
note that, due to Eq. (|2.51|) . 



/ Po[dtQ + C l , k u k ]dV = ^ n poQdvY 



and, assuming that the matter fields cancel on the boundaries <9D a , we get 



with 



^ [ J [Po(u\ + Q + OiuUdV) = y o f[dV, 



f[ = - p ip,i + Po Ck,i Uq- 

With the help of (|2.55|) and ([3.6)1 . this may be rewritten as 

Ml A 1 + F = 3 a + K a + 0(c" 2 ), 

where 



r = (i m ) 



po + po ( y - 3$ + n 



u l o + Pod > dV, 



(3.14) 



(3.15) 



(3.16) 



(3.17) 



(3.18) 



and 



jai= / (_ ai $._ po ^) d ^ 



'D 



(3.19) 



^3.20) 



Owing to Eq. ()3.7|) . Eq. ()3.17|) allows to compute the 1PN correction to the acceleration of 



the 1PN mass centers am: 



a (1) -a = x la c = — + 0(c ). 



(3.21) 



This is is the general form of the 1PN equation of motion for a perfect-fluid system in GR (in 
the harmonic gauge), according to the asymptotic scheme. To use this equation in practice, 
we must bring the integrals ()3.18|) - ()3.20|) to a tractable form, using relevant simplifications. 
This is done in the next Sections. However, we can already see that these integrals all 
depend on the Newtonian and 1PN matter fields, hence on the internal structure of the 
gravitating bodies. It is hence a priori clear that, unless a "miraculous" cancellation would 
occur, in (j3.21jl . of the different ways in which the internal structure influences the integrals 
(I3.18|) - ([3.20|) . the dependence on the internal structure should subsist in the final equations 
of motion. 
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IV. ACCOUNTING FOR THE GOOD SEPARATION BETWEEN CELESTIAL 
BODIES 



The "good separation" between the bodies of the system of interest means that the 



separation parameter 



?7o = max(r 6 /|a-b|) (r 6 = ~Sup xygD Jx - y|) (4.1) 

is small. We assume that the system of interest is described accurately enough by the 
equations of the asymptotic 1PN approximation, presented in the foregoing Sections — 
thus, a 1PN system S' is substituted for the "exact" system S. Then, to take the good 
separation into account in a clean asymptotic framework, we introduce again a (concep- 
tual) family of systems (S*), each being this time a 1PN system, and whh S< = S*. 
The family (S /T? ) is defined by initial conditions [12|. We refer to Ref. [12| for a more 
complete motivation of this approach, involving references to the literature. As pointed 
out there Q], this app roach seems to be new and came from the realization that, in 



our previous work 



11 1, the lack of considering the small parameter 77 within a such 



definite asymptotic framework (based on a conceptual family of well-separated systems) 
had led to the inappropriate neglect of some numerically significant terms in the 1PN 
equations of motion of our alternative scalar theory. Those terms turned out to be of 
order rf Q]. This order should be enough in the main solar system, for which we have 
7/0 = (radius of the Sun) /(minimum Mercury-Sun distance) ~ 1.4 x 10~ 2 , and 7/q ~ 3x 10~ 6 . 
In particular, using equations of motion derived from our alternative scalar theory and 
taking into account terms up to and including 77 s , we could reproduce the solar-system 
ephemeris DE403 of the JPL jsi| up to a 3"/cy difference (32I, the latter being due to the 
difference in the theories and in the approximation schemes. 



Let us state the initial conditions for system S' v . Owing to the 1PN equations, it turns 
out to be sufficient to define the initial zero-order density and velocity fields. To define 
Po(t = 0), we first define the initial position of the mass centers in system S" 7 : 

ag(* = 0) =ao(t = 0)770/77. (4.2) 

Equation ()4.2|) ensures that, at least near t = 0, the separation distances between bodies 
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are of order rj : 

(r° ab y = |ag - K\ = ordir,- 1 ) for a ? b. (4.3) 

Then, we just have to define the initial shape and size of each body (a) in system S" 7 , as some 
deformation of the initial shape and size of (a) in the system of interest corresponding to t)q. 
The simplest is to assume that in fact the bodies themselves do not depend on the separation 
parameter rj. This is expressed by the following definition of the density p^it = 0): 

pg(x,t = 0) =p (a + y,t = 0) if x = a[J + y with a + y G D a . (4.4) 

Equation (|4.4jl defines the field pg(i = 0) so that it is independent of 77 [setting 
Pg(x, t = 0) = if x does not have the form above for some a = 1, N]. 

To define the velocity Uq(£ = 0), we use the assumption that each body undergoes a rigid 
motion at the Newtonian approximation: 

u*=aj + fi£V'-fli), (n$> + Q$ = 0), or u = a + u, a A (x-a ), for x G D a . (4.5) 

This assumption is discussed in Appendix O It is shown there that this assumption is 
consistent to the accuracy aimed at in this work, i.e., to get the EMMC's up to the order 
rf included. We define the initial translation velocities of system S' v as 

(aint = 0) = ( V / Vo ) 1/2 ai(t = 0) (4.6) 

and the initial spin velocities by 

(^ ) nt = 0) = (i 1 /r ]0 ) 1/2 ^\t = 0). (4.7) 

Thus, the fields p^(t = 0) and u^(t = 0) are well-defined, and it follows from the 1PN 
equations that in fact all 1PN fields are then defined at t — 0. We assume that, as a 
consequence of this initial condition, the ord(?7~ 1 ) separation f)4.3|) holds true at any time in 
a relevant interval, in which one has moreover, consistently with ()4.6|) and ()4.7|) . 

(4)" = ord(r/ 1/2 ) (4.8) 

and 

(nW)' = ord(^). (4.9) 
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Assumption ()4.8j) for the magnitude of the Newtonian velocities is justified by the Newtonian 
estimate in a system with a dominating body, say body (N): 

$ « 2GM°/4 (4.10) 

together with the good separation ([4.3)1 . In the solar system, the angular velocities of 
the main bodies are quite small: at most of the same magnitude, in linear values, as the 
translation velocities — which is consistent with ()4.8|) and (|4.9|) . accounting for the fact that 
the size of the bodies is ord(?7 ). But in addition, the spins, including their axes, are 
very nearly constant. Since the spin evolution is determined by the rotational equations 
of motion, we cannot simply assume that the rates tl^ are very small, but have to derive 
it from these equations. In the asymptotic scheme used here, the zero-order (Newtonian) 
equations apply exactly to the zero-order fields, hence the the rates tl^) are determined by 
the Newtonian rotational equations. We shall prove in Appendix [B] that these rates are 
0(rj 3 ), provided that: i) the system is made of rigidly-rotating well-separated bodies, in 
the sense of Eqs. (|4.3|) . (|4.5jl . f)4.8|) and ()4.9|) . and ii) these bodies have a "quasi-spherical" 
inertia tensor, in the sense that 

= 0( V 2 ) (a = l,...,N; i,k = l, 2, 3), (4.11) 

where the 7 J - a ' ) 's are the eigenvalues of the inertia tensor 1^ defined in Eq. (|A7|) . This 
amounts approximately (or exactly, for homogeneous bodies) to assuming that the difference 
between the linear dimensions of each body in the different directions is 0(rj). It would be 
easy to modify the definition of the initial density field ()4.4|) so that it becomes compatible 
with (j4.11|) . e.g. by introducing an orthogonal affinity centered at a and leaving one 
dimension of the body unchanged while adjusting the other two. Thus, to account in an 
asymptotic framework for the numerical situation in the solar system, we need that our 
conceptual family of 1PN systems couples the good separation with the quasisphericity. 



(a) (a) 

7) -Tit 



In Appendix |X[ we use the orders in r\ given by ()4.3|) . ()4.8|) . ()4.9|) to evaluate the integrals 
(|3.18|) - (|3.2(J|) up to and including ord(?7 3 ) terms. [Of course, the density fields p , pi, etc., as 
well as the masses M° and M\ and other energy integrals, are ord(r] ).] There, we use some 



of the calculations done in Refs. 

n 



ii 



and also we use Newtonian or purely mathematical 
calculations of Fock [13j. We shall use the quasi-sphericity assumption (|4.11|) . leading to the 
estimate ()B8|) . only in two cases: to get Eq. ()A16|) giving I a \ and to get Eq. ()A44|) . 
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V. EXPLICIT 1PN EQUATIONS OF MOTION FOR WELL-SEPARATED, 
RIGIDLY-ROTATING, QUASI-SPHERICAL BODIES 



The explicit form of the 1PN correction to the equations of motion is got by inserting the 
explicit form of the integrals P, J a , K a , Eqs. (lAToT) . (1A351) . (jAll]l . (1A481) . and (lA32|) . into 
the general form (j3.21|) of the 1PN correction. In a first step, we obtain thus 



Xla 
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6T a + 2e a 
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b^a 

i 



x 16 -x la + 3 [(xi a -Xi 6 ).n° 6 ] n° 
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+ ^ab 

M(°).ao 
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T a 1 T a 

as 1 • a 1 -^lns "r ^21 
Jns ' 



£ab = 



T < +T 



+ 
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+ 0( V 7 / 2 ) + 0[c 



-2\ 



(5.1) 
(5.2) 



in which the space tensor Nl^ a \ related to the spin fiw, is given by Eq. ()A18|) . and the 
(spatial) vectors 1^, j£, L? ns and L£ ns are given by Eqs. (|ATT|) . (|Al?T|) . (|Al5|) and (|A47jl . 
respectively. The four latter quantities reduce to zero if the Newtonian density fields are 
assumed spherical in the sense of Eq. f|A20jl . It may be worth to repeat here that, in 
the asymptotic scheme which has been used in this work, the equations for the zero-order 



14. 



15. 



241 ] (see Sect. |TTJ). For this reason, 



quantities and the 1PN corrections are separated 
it is not directly possible to compare Eq. ()5.1|) with the LD equation. However, if we define 
the vector radius in terms of the full 1PN positions ()3.8|) : 



fab — |X a Xb| , 



n a b 



x a - Xb 

Tab 



(5.3) 



instead of defining it in terms of the zero-order positions, as in (jA36|) . then it is easy to 
derive the following 1PN expansion: 



n a b 



ab 



11 



+ 



1 f xi a - xib - 3 [(x la -xxb).n° 6 ] n° ab 



ab 



O (c- 4 ) . 



(5.4) 



We use this and the fact that, up to 0(rj 4 ) terms, one may use the spherical estimate (jA19|) 
for the external Newtonian potentials $^ a ^(a ) and ^-'(bo) in Eq. ()5.1|) . [This is independent 



27 



of the "optional" sphericity assumption (jA20j) . Moreover, we use the expression ()A32|) for 
M 1 in (j5.1|) and we switch to the notation ()3.8|) - (j3.9|) .] This enables us to group Eq. (j5.1|) 
with the Newtonian equation of motion (|3.1U|) so as to get 
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ab 
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ab 



to T a I j a 

ns _|_ i a _|_ 1 ns 2 1 

' Jns 



M° 



MP 



MP 



+ 0(r / 7 / 2 ) + 0(c- 4 ), 



(5.5) 



with, as before, matrix M^ a ^ of Eq. (jA18|) . and with 



d^b 



r bd (t = 0) 2 



3M 6 ° 



< 1. 



(5.6) 



The 5 b 's are constant up to 0(rf/ 2 ), see Eq. f|A13|) . Obviously, the term in (|5.5|) that 
involves this small constant may be absorbed, up to 0(c -4 ) corrections, in a redefinition of 
the zero-order masses, replacing them in ()5.5j) by 



M' h = M^ h {l + 5 b ). 



(5.7) 



Moreover, we recall that I° s , j° s , L? ns and L^ ns are given by Eqs. (|A17|t . (IA37j) . (IA45jl 
and ()A47)1 . respectively, and that these terms cancel in the "spherical" case in the sense 
of Eq. ()A20)1 . In that "spherical" case, the r.h.s. of ()5.5j) coincides exactly with the LD 
acceleration A^ D , apart from a term depending on the spin and on the internal structure: 



x„ - A 
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6T a MW 
K + Mf 
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ab 



MP ~mY) '~c? ] 



la 



c 2 M° 



[6(^ (a) ) 2 + nW.J]W] .x 0a + 0(r/ 4 ) + 0( 



(5.8) 



Here, 7 a is the spherical inertia moment (|A21j) . which does depend on the internal structure, 
and (uj^) 2 = is the square of the angular velocity. {O^.Q^ i s the "tensor" 

[see the footnote with Eq. (jA18|) ] with ik component fi^fi^ .} We note that this new 
term is order 3 in rj. It is clearly different from the two spin-orbit terms and the spin-spin 
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term appearing in the monopole-dipole equations of DSX |8j, e.g. because, on the basis 
of Eqs. (|OJ), (@~EJ) and (gU) here, the latter terms {Eqs. (6.32,33,34) in Ref. Q} are of 
order 7/ 4 , if and r] 5 , respectively. Many more structure-dependent parameters appear in the 
non-spherical case, see Eqs. flA37|) . (jA45|) and (jA47|) . But, as discussed in the Introduction, 
the departure from sphericity is unlikely to give numerically important 1PN corrections. 

VI. SUMMARY AND CONCLUSION 

We started from the family of initial data proposed by Futamase & Schutz Q|, which 
defines a family (S e ) of perfectly-fluid gravitating systems. This family is obtained by 
applying a Newtonian exact similarity transformation to the initial data defining the matter 
fields for the system of interest, S = S eo , and by assuming for the spatial metric an initial 
data which enforces spatial isotropy. That family allows to define the PN approximation 
in a clean asymptotic framework [14]. (A similar family, differing only in the initial data 
for the gravitational field, allows to do the same for an alternative scalar theory [3].) We 
call the corresponding PN scheme "the asymptotic scheme." By using a change of units 
depending on the parameter e Q], we are naturally led to postulate definite expansions 
for the independent fields, Eq. (|2.9|) for the matter fields and Eq. (j2.21|) for the metric. 
The latter is equivalent to the expansion of Weinberg j^j], but Weinberg did not expand 
the matter fields. However, in a first step, he did expand the energy-momentum tensor 
T, by assuming, for its different components, expansions in terms of the typical velocity 
in the system, his Eqs. (9.1.42-44). This is a slightly different expansion from the one 
we postulated in terms of the unique small parameter e, Eq. ()2.26|) here, but, together 
with the expansion of the metric, it leads to the same expanded equations: the Einstein 
equations, written in terms of the 1PN gravitational potentials, lead to Eqs. (|2.28|) - (j2.3U|) . 
while the dynamical equation leads to Eqs. ()2.31|) - (j2.34j) for tensor T. But since we do 
expand the matter fields, the explicit dynamical equations, as written in terms of the 
matter fields, are definitely different from those of Weinberg in his Sect. 9.8, which 
coincide with those obtained by Chandrasekhar Q|. The hydrodynamical equation for 
the 1PN corrections according to the asymptotic scheme valid for GR in the harmonic 
gauge was written here explicitly for the first time, Eq. (|2.59|) : Futamase and Schutz 
Eq. (4.28) 2 with (4.27) 2 -3 }, as well as Rendall {[2J], Eq. (44)}, did not write this equation 
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in a such fully-explicit form. Moreover, in these two works, equations for the "~PN" order 
were considered necessary. In the present work, we used e-dependent units in which the 
small parameter becomes e = c _1 and all fields are order zero (thus justifying the formal 
expansions in powers of c _1 ) and in which, moreover, only e 2 = c~ 2 plays a role in most 
equations — thus suggesting to postulate Taylor expansions in c~ 2 , for which there is no 
|PN term. In fixed units, this leads to an expansion of the metric and a general expansion 
of the energy- momentum tensor that are equivalent to those postulated by Weinberg jsoj l , 
involving merely even powers of e or merely odd powers. The calculations of his Sects. 9.1 
and 9.3 show that these expansions are consistent. 



We used the asymptotic PN scheme to derive 1PN equations of motion for the mass 
centers (EMMC's) in a weakly-self-gravitating system, according to GR (in the harmonic 
gauge). This was not done (with the asymptotic scheme) before. Like Fock Q, wm Q, 
and Brumberg [28], we define the mass centers through the rest-mass density, Eq. fQJ , 
and we obtain the EMMC's by integrating the 1PN dynamical equations in the domain 
occupied by the relevant body. One difference is that we use the asymptotic PN scheme 
instead of the "standard" PN scheme proposed by Fock and by Chandrasekhar [16J. 
As a consequence, we get separated EMMC's for the zeroth (Newtonian) order and for the 
first PN correction. In a first step, these are general equations valid for any 1PN system, 
Eqs. (I3.10j) and (j3.21j) . These general equations show that one has a priori to expect an 
influence of the structure of the gravitating body (e.g. the density profile) and its internal 
motion (e.g. the rotation velocity), since they depend on integrals of all matter fields. 
Another important difference with previous works (not only with Refs. Q, H, and Q), 
is that we use a definite asymptotic framework for the separation parameter rj, see Sect. IIVI 
which enables us to obtain results including an asymptotic error estimate. In the present 
work, we got explicitly the 1PN correction to the EMMC's when including the terms up 
to the order if included, Eq. (15.1)1 . To do that, we assumed, merely when calculating the 
1PN correction, that, at the Newtonian approximation, the bodies are rigidly rotating (as 
is nearly the case for the main bodies of the solar system), and quasi-spherical in the sense 
of Eq. ()4.11|) (which simulates in an asymptotic framework the real situation for these same 
bodies). 
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It turns out that the separated EMMC's of the zeroth and first order can be grouped 
together to give Eq. (|5.5|) . The latter may be readily compared with the Lorentz-Droste 
(Einstein-Infeld-Hoffmann) equation. In addition to new (structure-dependent) 1PN terms 
accounting for the departure from exact spherical symmetry, there is a new 1PN term 
depending on the internal structure and on the rotation velocity of the body considered, 
Eq. ([5.8)1 . It is worth to emphasize that the present derivation depends essentially on the 
weak-field assumption and is invalid for a strong-field system [although the definition of the 
mass center through the rest-mass density, Eqs. f)2.49|) and ()3.1|) . makes sense also for a 
strong field]. Thus, for a binary system of two neutron stars, each maximally spinning, the 
new term might lead to 10 % corrections to Newtonian physics, whereas in fact, if the binary 
is widely separated, the Newtonian dynamics is known to be a good approximation. 
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APPENDIX A: THE INTEGRALS P, J a AND K a FOR WELL-SEPARATED, 
RIGIDLY-ROTATING BODIES 

Let us first note the following correspondence between the notations of Weinberg j^ol ]. 
used in the present work, and those of Fock, also used in Refs. jlO, 11, l^J {see Eqs. (|2.3U|) . 
(l2~3oT) and (l2~oDl) here, and Eq. (73.03) in Fock Q}: 



-U, Q 



-4Ui. 



(Al) 
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1. Integral I a 



This integral [Eq. 1)3. 18|) ] is equal to 



po + po [ y - $ + n 



u l dV 



-2 [ {- P ^v i Q )dV + 2 f (- Po( f )a ui ) )dV+ [ poQdV 

JD a JD a JD a 



= I™ + 21 f + 21 f + I ai 



1 ' 



(A2) 



.0 



Integral If on the r.h.s. was in fact denoted by I m in Refs. |lOUlll.ll2j] [it indeed played the 
same role there as the integral (I3.18J1 plays in Eq. (|3.17|) here] and it is given in two pieces 
by Eqs. (A20)-(A2lV in Ref. the second piece, Eq. (A21) in Ref. [111 ], is just integral 
If above. In Ref. JjJ , it could not be assigned a definite order in i] to the remainders, but 
now it is in general easy to do that, using Eqs. (|4.3jl . ()4.8j) and (|4.9|) . 5 We have thus: 



If = If + (M a °a*/2 + 2T a + 4e a )d* + M + + + «(r/ 7/2 ), (A4) 
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'D a 



(A7) 
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- 1 - a 



a .) 



Po^aix 1 - a J )dV. 



Da 



(A8) 



5 An exception is for the order in Eq. 1A4|) here. It was proved in Ref. [^3, Eq. (A12), that Fock's 
"Lyapunov's equation" {0|, Eq. (73.26)} indeed applies to a body in a well-separated system, with a 
definite remainder in 77: 

PO+P O n O -Po(^a-0a) + O(?7 3 ), (A3) 

where f2 a is defined in Eq. (| AT|) below. [In this proof, the fact was used that 0$ = 0(i] 3 ) as a consequence 
of the Newtonian rotational equations — this is proved in Appendix [B] below, yet that proof needs quasi- 
spherical bodies in the sense of Eq. g~TTJ.] We have = ord^ 1 / 2 ) [Eqs. l|P |l -l|P fl ]. Equation l|Al)l . 
including the order of the remainder, follows by a straightforward computation. 



32 



The third integral is immediately computable from Eq. (|4.5|) above, in terms of the potential- 
energy integrals (jA6|) : 

If = [ (-aA4) W = 2(e a aj + e^?). (A9) 

By dividing the potential Q (the negative Newtonian potential associated with the density 
A T t0 = 4p Mo) i 11 ^ "external" and "self parts and ( a i, exactly as for $ in Eq. (|3.11j) . 
one finds [cf. Fock's integrals (75.36) and (76.29-30)]: 

If = -41 f + If, (A10) 

with 



If = [ Po zl a) dV = M a Zi a \z) + Otf/*) 
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Summing these different contributions, we get 
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(M^/2 + 2T a )d* + (aSfiff/W + 2T aj )^? 
-3M a V $ (a) (ao)-4GM Q °E Mb °° 



|a — b 



o 



-3/^^(ao) + 4^ (A12) 

fe _, a Ca |a -b | 



When differentiating this with respect to the time, as is needed for insertion in the equation 
of motion (|3.17p . we note that, for well-separated, rigidly-rotating and quasi-spherical bodies 
(at the Newtonian approximation), the rate of the (Newtonian) angular rotation velocity 
is 0(i] 3 ) [Eq. (jB8|) ]. Moreover, independently of the quasi-sphericity assumption (J4.11j) . we 
have by the Newtonian rotational equations (jB3|) {cf. Fock [13], Eq. (72.32)} 

dT a 



dt 



=0{rf/% (A13) 



In addition, due to the assumed rigid zero-order motion (|4.5|) . the (Newtonian) inertia tensor 
is constant in the frame following the (Newtonian) rotation of the body so that its time 
derivative in the starting harmonic coordinates is known explicitly as 



>ody, 





i%=^I$+ti$^- (M4) 
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In the same way, we have by (jB8|) and Fock's Eqs. (72.24) and (74.06) 

<1T^ _ 1 (a) (a) (a) (a) {a) (a) (a) (a) , r(a) f , i i\( T k k\( I I 

J Do, 



dt 

Therefore, we get 



d„ 

(A15) 



F = M°(^.ao)^ + {M°[^ 

-AGM®^- Y\ - M(a) -ao + K s + 0(^ 2 ), (A16) 

dt j—^ a — b 

b=£a 

C - ^/jf^^ (A17) 

f>7^a I | 

in which the rates of the inertia tensor are given by (jA14|) . and where 



(i 



MW-ao = (M^d l ), - nS } /if Off- (A18) 

To have F up to the 0(r] 7 ^ 2 ) remainder, we may use in (jA16j) and (|A17|) the spherical (or 
monopole) estimates of the Newtonian potentials; in fact the dipoles cancel due to Eq. (|H.4|h . 
so that 

$(a) ( a o) = -E^% + °(^)' ^ $(a) (-o)] = 4E^b + °^ /2 )- (^19) 
r~ a o — D o at TT a o ~~ D o 

In the case that the Newtonian mass density p is assumed spherical inside each body: 

Vx € D«, p (x) = pair), r = |x-ao| (a = 1, N), (A20) 

we have 

/ \ 47r r Ta 

I^=laS lk , 7a = y/ n r*p a (r)dr. (A2I) 

Q, 



Then the rates (|A14|) cancel. The same is true for the T a j's jl3j, hence I™ cancels. 



6 Recall that we are using a fixed system of harmonic coordinates x* 1 , such that the initial condition has the 
form Q2.1 |l - (|2.4(l . However, the equations written here are covariant under a "Cartesian" spatial coordinate 
change, 11 = RijX^ with R = (Rij) a (constant) orthogonal matrix, R € 0(3, R). Hence the 

upper or lower indices need not be distinguished (they may be exchanged by the Euclidean metric that 
has components Sij in all these coordinate systems). Moreover, these equations may be brought to a 
form covariant under an arbitrary spatial coordinate change: y° = x°, y 1 = ijf{x 3 ), provided some care is 
taken, e.g., in handling the time derivatives: the components of the velocity vector v along a trajectory 
t i ► x(i) are, in any space coordinates y z , the time-derivatives v % = y 1 , but the acceleration A must be 
defined as the "absolute" time derivative of the velocity vector using the connection associated with the 
Euclidean metric; the components of A are x! % only in "Cartesian" coordinates x n . 
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2. Integral J 



To calculate integral J ai [Eq. (J3.19|) ]. we introduce Fock's auxiliary potential 

W(x, t) = j G |x - y| p(y, t) dV{y)/2 (A22) 

and we note that, due to Eqs. ()2.29|) and (|2.56|) . and since AW = — $, the potential if) is 
given by 

-if, = B + d 2 W/dt 2 } (A23) 
where B satisfies the Poisson equation 

AB = -AnGat. (A24) 

Therefore, integral J a% is exactly the integral denoted so in Ref. ^(|, Eq. (4.11) and after, 
though here the field en is defined by Eq. (12.561) . This integral has been studied in detail 

nn 

in Refs. |lfl | -[12 j . under the assumption that the boundary condition for B ensures that is 
is indeed the (positive) Newtonian potential associated with the matter field a±, the latter 



having spatially-compact support. We have {[12], Eq. (A14)} 

Ox^ x=ao 



jai _ L a 4 = GM^^-r V 



h Pbj; r^I 



- — bo I |x — b 

- a 'lH- p "HH +0 ^' (A25) 

with 

r d 3 W 

a a = [ aidV, Paj= f *i(x)(x> -a&)dV(x). (A27) 
The coefficients a a and /3 a j (a = 1, N, j = 1, 2, 3) defined by (|A27|) have been computed 



in Ref. |12|, though with a slightly different o\ from the one valid for GR (in the harmonic 
gauge), Eq. (|2.56|) here. The calculation is easy, using the rigid velocity field (|4.5|) together 
with Eq. (jA3|) and a Taylor expansion of at ao- Due to the presence in ()2.56|) of the 
field po, we have now to know the following integrals [Fock's Eqs. (74.24-25)], which also 
follow from "Lyapunov's equation" (|A3|) . and which have the same remainder as the latter: 



3 / Po dV = e a - 2T a + 0(r, 6 ), (A28) 
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where 



2 / Po (x j - og) dV = rjai ~ T aj + 0( V 3 ), Vaj = (B« + i^)/2, 



We find thus: 
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and 
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9T a + ^-e a + 0(r] 3 ) 



(A33) 



/3 aj = M^ai - ai) + tariff + 3T aj + Vaj + 0( V 2 ) = M a Vi - a ) + Vaj + 0(rj). (A34) 



To the purpose of getting the EMMC's up to rf included, the second estimate of j3 a j, up to 



0{rf) not included, is enough since f3 a j is multiplied by 0{rf) in ()A25|) . To the same purpose, 
we may use the spherical estimate (jA19|) i for $( a ) in the expression (jA33|) for a a , and also 
in the expression f)3. lUj) of the Newtonian acceleration a , which enters the PN correction 
to the acceleration (j3.21|) . We also use Eqs. ()3.7j) and ()3.8|) to reexpress the terms with the 
P 's. This gives 



fa j a 



M^ao 



b^a 



a b + (a a - Ml 



MS 



n 



ab 



ab) 



+ \ a 
1 Jns 



G^t^ [x 16 -x la + 3 ((x la -x 16 ).ny n° ab ] + 0( V 4 ). (A35) 

b^a ^ Tab ' 



7 In this integral, the integrand, say /, is not Lebesgue-integrable, but we have / = / sp her + Sf, where 
/sphcr is the integrand corresponding to spherical symmetry (with </> a (x) — <p a (r), r = |x — ao|, and 
4>a{r) = —GM®/r for r > r a ), and where Sf is Lebesgue-integrable. Hence, the integral exists in the sense 
of Xpacc / dV = lim «-oo J| x - ao |<fl / W, because J| x _ ao |< B /sphcr dV = for any i? > 0. 
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Here, 



- ah = |a - b 1 = |x 0a - xo6| , n a6 = g (A36) 

r ab 



and 



J ns 



The latter cancels, as do the rjbj 's, when all bodies are spherical in the sense of (jA20|) — see 
Eqs. (JMH|) 2 and (IA42J) . 



Integral ()A26|) was calculated in Ref. 11( to an order which turns out to be sufficient. It 
is given by {[ll|, Eqs. (A14), (A20)}: 



l " = At jj'°3¥ik AV + / D /o^dV + (c-^Lf + ^ + 0(c-), (A38) 

where the "self part is given by the following equation {[11], Eq. (A21)}, which is an exact 
one: 



Lf = Tt (-eA + B^'dS - e^df + flWfiW) , (A39) 

Q 

B$ = I [5 lk {V<p a ) 2 /2 - <f> a ^ k } dV/AnG. (A40) 

J space 



Due to the assumed rigid motion (|4.5|) . the rates of the quantities £ a , e a j, Bj) k are 
known: e a = 0, follows the rule (|A14|) . and we have [k| 

e aj =^e akj B% = ^B^ + ^Bfi + ^B^. (A41) 

Moreover, in the case that the Newtonian density po is exactly spherical in body (a), then 
we have |ll[ 

e aj = 0, B^ = 5 lk e a /3, B% = 0. (A42) 

Recall that we shall obtain the EMMC's for a family of well-separated and quasi-spherical 
bodies, i.e. a family of bodies, indexed by the separation parameter 77, and whose family all 
bodies become closer and closer to being spherical as rj — ► 0. For a such family, we should 
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have 



and Bj) k — > as i? 



0. 



Therefore, we get from ()A39|) and ()B8|) : 



with 



We have L"* ns = (exactly) if the Newtonian density fields are spherical in the sense of 
Eq. (lA~20l) . 



l(a) 



□ (") n( a ) 



(A43) 



(A44) 



(A45) 



The "external" part is given by Eq. (A24) of Ref. [llj, of which the last term is 0(t] 4 ), 9 
hence may be omitted here: 

b ol ijkyjd 3 l a o ~~ b oi 



Lai 
2 



b^a 



da l daQ 



J 2ns 



2 ^ 3 

b^a 



da^daQdaQ 



d 3 |a - b 1 

dalda^dal 



+ Lf m + 0( V « 



(the latter cancels in the case of spherical symmetry), or more explicitly 

3«,b ) 2 - hi 



gm; 



-J2 M i 



b^a 




n L b o n° 6 - b c 



ab 



+ 



n° a6 -2(n° a6 .b )b Q 

ab) 



(A46) 



(A47) 



(A48) 



3. Integral K a 

To calculate integral K ai [Eq. (|3.2(J|) ]. we first note that, due to the Poisson equation 
(JPty with the expansion fl23SJi, we have {Eq. (75.27) of Fock |l3|}: 

/ PoCka^u'dV = 0. (A49) 



This does not strictly follow from the assumption Q4.11[l . which concerns merely the departure from 
sphericity of the inertia tensor. However, if one would explicitly define, in a natural way, the density 
field p[j satisfying (|4.11|) . as outlined after that equation, then he should indeed obtain (|A43(1 . and even 
probably the e a j 's and the Brf k 's should be 0(rj) or higher. To be logically consistent, wc may content 
ourselves by simply assuming that the family does satisfy l|A43|l in addition to l|4.11JI , 
because the assumed rigid motion implies that = 0(ftij) — O^ 1 / 2 ), see Eq. (|A14p . 
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Hence it merely remains the "external" part. I.e., we have 



K° 



P oZt}u k dV. 



(A50) 



To compute it, we use Fock's Eq. (76.23), which we rewrite in the present notation, also 
giving the (easily-evaluated) order of the remainder: 

d 1 



^ Q) (x) 



. 4GE MX 



X - bn 



4G 



lj dxi |x - b n 



+ 0(r] 7/2 ). (A51) 



bj^a b^a 

Inserting ()A51|) and the Taylor expansion of 1/ |x — b | into Eq. (jA50|) with the rigid velocity 
field l|P|) . we get 

K a = AGM° a £ M°a .b -% + 0(r? 4 ). (A52) 

b+a ^ 

APPENDIX B: NEWTONIAN SPIN EVOLUTION FOR A SYSTEM OF WELL- 
SEPARATED, RIGIDLY- ROTATING, QUASI-SPHERICAL BODIES 



We start from Fock's Eqs. (72.06) and (72.09), which, combined together, write 



Po 



0®^ , k 

(l'n)-^-r- + {X K 



dx k 



dx i 



dV, 



where 



Po [{x l - a )u k - {x k - a k )u ] dV. 



(Bl) 



(B2) 



Equation (jBl|) may be seen as the application to body (a) of the Newtonian theorem stating 
that the rate of change of the total angular momentum of a system is equal to the sum of 
the external torques on the system. Recall that the Newtonian equations apply exactly to 
the zero-order quantities (Subsect. Ill Cj) . and that we are actually considering a family (S /?? ) 
of 1PN gravitating systems (Sect. HVjl (although we omit the superscript r] on the fields for 
icitv 

3, 

Eq. UOJ here}: 
d 



the simpli city of notation). Due to the good separation, the r.h.s. of (jBl|) may be evaluated 
as {Fock 13j, Eq. (72.13) in which we evaluate the order of the remainder with the help of 



E 

b^a 



3GM 6 Vo - K) 



ao 



b k )I^ - U - b n )& ] 



(B3) 



As to the angular momentum tensor (jB2|) . using the rigid velocity field ()4.5|) as well as the 
definitions of the Newtonian mass center (j3.4|b and the inertia tensor (jA7|) i . it is easily 
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calculated to be 

M^=n$I^-^I$. (B4) 

We differentiate this with respect to the time, using the rate (jA14j) of the inertia tensor. 
Then, we rewrite the result in space coordinates (x n ) which are Cartesian for the Euclidean 
space metric (that one having components 5ij in the given harmonic coordinate system (x M ) 
utilized) and which, at the current time t, bring the inertia tensor of body (a) to the diagonal 
form (the x h 's are deduced from the x l 's by a space rotation): 

4° } = Vrj'V (B5) 

This yields 

^ = ng) ( T f + 7 w) + nwnw ( 7 w _ 7 w) . (B6) 

We also simplify (jB3|) with the help of (jB5|) . and we equate the result with (jB6|) . We thus 
get in the coordinates (x'*): 



(a) _ (a) 
rS(a) _ 7j ^fc 
ik ~ (a) , (a) 

Ti +% 



v- 3GM°(a -6 )(ag-6 fc ) _ (a) (a) 
l a o ~~ D o| 



+ Otf). (B7) 



If, now, we account for the fact that = 0(r// 2 ) [Eq. 1)4.9 J) ] and for the quasi-sphericity 
assumption (J4.ll)) . we do see that 



Oi? = 0(^ 3 )- (B8) 

APPENDIX C: JUSTIFICATION OF ASSUMING A RIGID ROTATION FOR 
WELL-SEPARATED BODIES 

Let us shortly discuss the possibility for a perfectly-fluid body in a weakly self-gravitating 
system to have a rigid motion at the zero-order approximation [Eq. ()4.5|) ]. Using the continu- 
ity equation 1)2.51)1 . we rewrite the zero-order dynamical equation ()2.52jl as Euler's equation, 
and in the latter we insert the rigid velocity field ()4.5j) . This yields 



a + (n« - n$?ng) & - 4)] = - Po $, - Po ,. (ci) 



Thus, the r.h.s. must depend linearly on the position x, as does the l.h.s.. According to 
Eq. (jH.lljl . the Newtonian potential <3? decomposes into the self-potential <ft a and the external 
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potential $(°). If the body is isolated, the latter cancels. In that case, Eq. flCl|) can certainly 
be verified, for it is well-known that an isolated rotating mass made of a perfect-fluid body 
is dynamically possible (may be under certain restrictions on the state equation). Thus, 

po (njf - ng'nS?) (* j - ^ = - Po <p a , t - Po>i (C2) 

can be satisfied exactly. (The mass center has no acceleration in that case, of course.) 
The presence of external bodies produces time-varying tidal forces which prevent the exact 
validity of the equilibrium (|C1|) . This means that in fact the body will undergo some time- 
dependent deformation: tides. If we consider a family of well-separated systems, we have 



^ a) (x) = ^^(ao) + $g(ao)(^' - 4) + 0(r? 4 ) (x 6 D a ). (C3) 

Inserting this Taylor expansion into the equation of motion for the zero-order mass center, 
Eq. (j3.1(Jj) . and accounting for (|3.4jl . we find that 

K = <S>f(z ) + O( V % (C4) 

We can then rewrite (jCl|) as 

Po (n<? - nS } ng) (*>' - oj) = -p + *«(ao)(^ - 4)] - p 0ji + Ofo 4 ). (C5) 

This shows that, to accommodate a given spin (fi^), we cannot start from a solution 
of Eq. (|C2|) and modify just the spin rate [which, we recall, is primarily subject to 
Eq. (|B7|) ]: we must also modify the pressure field po, which (through the state equation) 
determines the field po, the latter determining in turn the self-field <p a . This would need a 
detailed study. But if in the equilibrium equation (|C1J1 we neglect not only 0(rj A ) but even 
0(r] 3 ), then we are left with 

po (n« - nifng) - 4) = - POil + o( v 3 ), (ce) 

in which, owing to Eq. (IB 811 . we may further neglect Q ? — provided the system is quasi- 



spherical in the sense of Eq. (j4.11|) . This is all what we used. In particular, this is what 



was used in Ref. |33( to derive "Lyapunov's equation" reproduced above, Eq. ()A3|) . Thus, 
to get the translational equations of motion up to rf included, we need only to solve the 
Newtonian internal equilibrium up to the order if included. The latter is not influenced by 
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the external bodies, Eq. (|C6[) . hence it is compatible with a rigid rotation (at the Newtonian 
approximation), as postulated in Eq. ([4.5)1 . 
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